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We model the electrons on a monolayer graphene in terms of the compact and non-compact U(l) 
lattice gauge theories. The system is analyzed by the strong coupling expansion and is shown to 
be an insulator due to dynamical gap formation in/around the strong coupling limit. This is similar 
to the spontaneous chiral symmetry breaking in strong coupling gauge theories. The results from 
the compact and non-compact formulations are compared up to the next-to-leading order of the 
strong coupling expansion. Excitonic modes and their dispersion relations in the insulating phase 
are also investigated: it is found that there arises a pseudo-Nambu-Goldstone mode obeying the 
Gell-Mann-Oakes-Renner type formula. 
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1. Introduction 

After its first experimental discovery in 2004 [j|], graphene (monoatomic layer material of 
carbon atoms) has widely attracted theoretical and experimental attention [Qj. Due to its hexagonal 
lattice structure, charge carriers on a graphene reveal a linear dispersion relation around two "Dirac 
points" in momentum space [Q], so that quasiparticles at low energies can be described as two 
species of massless Dirac fermions in (2+l)-dimensions [Q]. The symmetry between two triangular 
sublattices of graphene is referred to as the "chiral symmetry". 

In the low-energy effective theory of graphene, the effective Coulomb interaction between 
charge carriers is enhanced by 300 times due to the small Fermi velocity v F . Such a strong Coulomb 
interaction may turn a suspended graphene in the vacuum from semimetal to insulator by the for- 
mation of a finite spectral gap of quasiparticles [0]. This mechanism is similar to the spontaneous 
chiral symmetry breaking and associated fermion mass generation in quantum chromodynamics 
(QCD). Various attempts have been made so far to study the gap formation in monolayer graphene 
by using the Schwinger-Dyson equation 0, g], the 1 /N expansion [[| [l(J pQ , the exact renor- 



malization group []12[], and the lattice Monte Carlo simulations \\13[ \14\ \l% |lq ]. These works are 
mainly focused on the critical region of semimetal-insulator transition or on the behavior for large 
number of flavors. 

In this work, we rather focus on the strong coupling region of the system and study the low- 
energy effective theory discretized on a square lattice jl8| , [l9[ l. By using the strong coupling expan- 
sion of the compact and non-compact formulations of the gauge field, we study the gap formation 
due to the spontaneous "chiral symmetry" breaking as well as the dispersion relations for collective 
excitations. Typical energy scale of the emergent excitations are also estimated with the use of the 
intrinsic length scale of the original honeycomb lattice. 



2. Low-energy effective theory of graphene 

2.1 Effective action in the continuum limit 

With the annihilation operators of the electrons on the two triangular sublattices of graphene 
(a a and b a ) near the two Dirac points K±, we can construct a 4-component spinor in the momentum 
space, V'ct(p) = (a CT (K + + p),& CT (K + + p),& f7 (K_ + p),a f7 (K_ +p)) r . Here a =|,| denotes the 
original spin of the electrons. The Euclidean effective action for graphene is then written as [f7], [T^] 

S E = £ / <h® % (D[A] + m) y a + £ / dx^ (djA 4 ) 2 , (2. 1) 

where the Dirac operator reads D[A] = 74(^4 + 1A4) + v F Y*v=\z Yi^i- This is analogous to the action 
in "reduced QED" [pQ], in which the fermion y in (2+l)-dimensions is interacting with the U(l) 
gauge field A in (3+l)-dimensions. The Hermitian / matrices obey the anticommutation relation 
{7^1 7v} = 25^v The gauge coupling constant is defined by g 2 = 2e 2 /£o(l + £)> with the electric 
charge e, the dielectric constant of vacuum So, and the dielectric constant of substrate e. Due to the 
small Fermi velocity of quasiparticles, v,, = 3.02 x 10~ 3 , one may adopt the "instantaneous approx- 
imation" in which the spatial component A is neglected. From the absence of the z-component in 
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the Dirac operator, this model possesses a continuous global U(4) symmetry generated by 16 gen- 
erators {1, 73, 75, 7375} <8> {1, a}, which is the extension of the continuous U(l) charge symmetry 
and the discrete Z2 sublattice exchange symmetry of the original honeycomb lattice. The explicit 
chiral symmetry breaking term is represented by the mass m in Eq.(2.1). 

After performing the scale transformation in the temporal direction, X4 — > x^jv F , A4 — > v F A^, 
the Dirac operator becomes independent of v F . This scale transformation changes the mass m into 
the effective mass m if = mjv F , while the Coulomb coupling strength is enhanced as g% = g /v F 
which is about 300 times larger than that of QED. Since the inverse effective coupling strength j8 = 
l/gl is 0.0369 in the vacuum-suspended graphene, the expansion by j3 (strong coupling expansion) 
would work well. 



2.2 Regularization on a square lattice 



We discretize the low energy action Eq.(2.1) on a square lattice with a lattice spacing a. Since 
the original honeycomb lattice has a lattice spacing a Hc ~ 1 .4A, we make an approximate iden- 
tification, a ~ a Hc , so that we can carry out the strong coupling expansion. The quasiparticles in 
monolayer graphene are described with a single staggered fermion %, because its eight doublers 
can be identified as 4(spinor components) x 2(spin) degrees of freedom. As a consequence, the 
lattice action for fermions on graphene is written as [ |l4| ] 



1 



I (V+( X )-V-(x))+m*M(x) 



M=l,2,4 



(2.2) 



withfermionicbilinearsM(x) = %{x)%{x),V+(x) = r\ tl {x)x{x)U yl (x)x(x+ ft),V^ (x) = r\ yL {x)x{x J r 
fi)Ul(x)x(x), where /I = 1,2,4. The staggered phase factors corresponding to the Dirac 7- 
matrices are rj 4 (x) = 1,tJi(jc) = (— l) A4 ,T] 2 (x) = (—l) X4+Xl . is the U(l) link variable, where the 
temporal link is defined as U^x) = exp [/0(x)] (—71 < 6 < n), while the spatial links C/i,2,3 are set 
to unity as a result of the instantaneous approximation. In the staggered fermion formulation, the 
global chiral symmetry U(4) shrinks to U(l)v X U(1)a, with the ordinary U(l)v charge symmetry, 
and the axial U(1)a symmetry generated by e(x) = ( — 1 ) A " 1 + ' V2 + - V4 . 

As for the pure gauge action, we consider two types of formulation. One is the compact 
formulation which consists of plaquettes made of U(l) compact link variables: 



4 C) =i8l I l-Re(u 4 (x)Ut(x + j 



*M 7=1,2,3 



(2.3) 



The other is the non-compact formulation with the gauge angle 6: 



,(NC) 



?E I [e(x)-e(x+?)] 2 . 



(2.4) 



' *M 7=1,2,3 



The compact formulation has photon self-interactions which are absent in the non-compact formu- 
lation and in the continuum theory. 
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3. Strong coupling expansion 



Expanding the partition function Z by the inverse coupling strength j3 and integrating out the 
link variables order by order, we obtain the effective action S x in terms of fermions [[P 



J[d X dx}[de] 



[d%d%\e 



(3.1) 



Since the link integration selects the terms in which the link variable cancel with its complex 
conjugate, various 4-fermi interaction terms are induced as shown in Fig.[jJ With the compact 

formulation, we obtain the leading order (LO) and the next-to-LO (NLO) effective action, Sy and 

e (i)c 



respectively, as 



,(i)c 



I 



o E [Vf(x)-Vf(x))+m*M(x) 

Z 7=1,2 



-£M(x)M(* + 4), 



EE V+ (x)Vj- (x + 4) + Vf (x)V+ (x + 4) 

*(3) j=l,2 



(3.2) 
(3.3) 



Since the pure gauge term So vanishes in the strong coupling limit (/3 = 0), the compact formulation 
and the non-compact one give the same result in the LO. In the NLO, the effective action from the 
non-compact formulation, sjjj NC , is twice that from the compact one, S^ C , so that observables of 
the both formulations are related as ^ NC (j8) = <^ c (2/3). 



M(x + 4) CCM 



Vf(x+4) (( 



M(x) (gQ 




(a) 



(b) 



Figure 1: Induced four-fermion interaction in the strong coupling expansion. The open and filled circles 
represent fermion fields % and %, respectively, (a) In the LO, the time-like links (red arrows) in the fermion 
action Sf cancel with each other to leave a spatially local interaction, (b) In the NLO, the time-link links in Sf 
are canceled by the time-like links in a plaquette Sg (blue arrows) to leave a spatially non-local interaction. 

In order to linearize the induced 4-fermi terms and to integrate out the fermions, we introduce 
complex bosonic auxiliary fields by the Stratonovich-Hubbard transformation. As for the LO 
4-fermi terms in Eq.(3.2), we introduce an auxiliary field (f>(x) = <t> a {x) + i&{x)fy n {x). Another 
auxiliary field A = X\ + ihi is introduced to linearize the NLO terms in Eq.(|3.3|); the mean field 



value of A is determined by requiring the stationary condition of the effective action. Then we 
arrive at the LO and the NLO effective potential (free energy) written in terms of as 



p(l)C 
eff 



4 ivi 2 
4 7=1,2 



InfG-HM)], 



G(k;0) sin 2 k } 



(3.4) 



(3.5) 
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with the bosonic effective propagator defined as G 1 (k;0) = \<j>/2 — m*| 2 + L/=i,2 sin 2 fcy and the 
momentum integration as / k = {2n)~ 2 j* n dk{ f 7C 7[ dk 2 . Fig.|| shows in the chiral limit (m = 
0): Its minimum corresponds to the chiral condensate (the order parameter of chiral symmetry 
breaking) a = \(XX)\ = K0)l> so tnat the spontaneous "chiral symmetry" breaking takes place in 
the strong coupling limit. The effective mass of fermions reads Mp = m + (v F /a)(aa 2 /2), from 
the mass term of the effective action. 



4x1 3 : 



o 

2x10 



0.0 




-2x10 



Figure 2: The free energy in the strong coupling limit, Fg (<j>), in the lattice unit as a function of | (j) |, in the 
chiral limit (m = 0). 

Since is an increasing function of |0 1, the chiral condensate o drops as /3 grows. In other 
words, the chiral symmetry tends to be restored as the coupling strength becomes weaker. Up to 
the linear terms in /3 and m, we can calculate a with the compact formulation as 



a c (j8,m) ~ (0.240 -0.297/3+ 0.0239 ma)a 



(3.6) 



As we mentioned, the condensate in the non-compact formulation is simply obtained as a NC (j8) = 
a c (2/3 ) up to NLO. The behavior that a NC (j8 ) drops faster than a c (j8 ) is consistent with the results 
of the Monte Carlo simulation for the same lattice model [16]. Taking aT^ ~ az} = 1.39 keV as 



a typical cutoff scale of our system, we obtain (7 c (/3,m) ~ [(0.680 - 0.421/3 + keV] . The 

dynamical fermion mass is estimated from the value of the chiral condensate as Mp ~ (0.523 — 
0.623/3) eV + 3.05m, which is much smaller than the momentum cutoff scale ~ 7iv F /a = 13eV 
of the original honeycomb lattice as long as the bare mass m is small enough. 



4. Collective Excitations 

Here we study the fluctuations of the order parameter (j>(x) around the symmetry broken state 
(0) = —a: the phase fluctuation ("7r-exciton") corresponds to ty n {x) and the amplitude fluctuation 
("(7-exciton") corresponds to <j> a (x). Propagators of those modes, K are derived from the second 
derivative of the effective action Self [<p] with respect to the corresponding fields ty a ,n- Their excita- 
tion energies are obtained from the imaginary pole of the propagator, D7 (p = 0,(0 = iM a>n /v F ) = 
0. 

As for the 7i-exciton, we obtain a mass formula in the leading order of m as 

2v K I m 

M * = vvW^oT <41) 
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Since M n vanishes in the chiral limit (m = 0), this mode serves as a pseudo-Nambu-Goldstone 
(NG) boson emerging from the spontaneous breaking of chiral symmetry. From the axial Ward- 
Takahashi identity corresponding to the infinitesimal local chiral transformation, a simple relation 
can be derived, (F^M n ) 2 = mo, which is analogous to the Gell-Mann-Oakes-Renner relation for 
the pion in QCD [gl]]. Here the temporal "pion decay constant" is defined by the matrix element, 
(0|/f ial (0) = 2Flp\, with the axial current Jf a \x) = |e(x) [V~ (x) - V+(x)] . 

By solving the pole equation numerically, the a-exciton is found to be a massive mode with 
M G ~ (1.30 — 0.47/3)(v f /a) + 22. 6m. Since M c acquires a large value comparable to the cutoff 
scale Ea, application of the low-energy effective theory in this channel is not quite justified. 



5. Conclusion 

We have investigated the behavior of monolayer graphene analytically in/around the limit of 
strong Coulomb coupling, by means of the strong coupling expansion of U(l) lattice gauge theory. 
As for the pure gauge action, we have compared the results from the compact and non-compact 
formulations. In either case, we find that "chiral symmetry" (the sublattice exchange symmetry 
in the original honeycomb lattice) is spontaneously broken in the strong coupling limit with equal 
magnitude of the chiral condensate. As the coupling strength becomes weaker, chiral condensate 
from the non-compact formulation drops faster than that from the compact one. These results up to 
NLO in the strong coupling expansion agree qualitatively with the extrapolation of the numerical 
results of the lattice Monte Carlo simulations. 

We have also examined the collective excitations associated with the chiral symmetry break- 
ing in our approach and have derived their dispersion relations. The phase fluctuation of the chiral 
condensate, the "7T-exciton", behaves as a pseudo-NG boson, like the pion in QCD. Experimental 
observation of such mode in vacuum-suspended graphene would be a good evidence for the spon- 
taneous chiral symmetry breaking in the strong coupling regime. The amplitude fluctuation of the 
chiral condensate, the "a-exciton," acquires a large mass comparable to the intrinsic cutoff scale 
E\, so that it needs further investigation without the low-energy approximation. 

There are several directions to be investigated in future: Behavior of the present model on a 
square lattice at finite temperature and finite chemical potential still remains an open problem both 
analytically and numerically. Formulating the strong coupling expansion on a honeycomb lattice 
would be of great importance. Extension of our strong coupling approach to the analysis of bilayer 



graphene, which has been attracting attentions recently [22], would be also of interest. 
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